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After surveying the experimental evidence for concentration coupling in the shear banding of 
wormlike micellar surfactant systems, we present flow phase diagrams spanned by shear stress (or 
strain-rate) and concentration, calculated within the two-fluid, non-local Johnson-Segalman (d-JS- 
(j)) model. We also give results for the macroscopic flow curves E(7, (^) for a range of (average) 
concentrations (j). For any concentration that is high enough to give shear banding, the flow curve 
shows the usual non-analytic kink at the onset of banding, followed by a coexistence "plateau" that 
slopes upwards, dT,/d-y > 0. As the concentration is reduced, the width of the coexistence regime 
diminishes and eventually terminates at a non-equilibrium critical point [Ec, 0c, 7c]- We outline 
the way in which the flow phase diagram can be reconstructed from a family of such flow curves, 
E(7, (j)), measured for several different values of (j). This reconstruction could be used to check new 
measurements of concentration differences between the coexisting bands. Our d-JS-(;!> model contains 
two different spatial gradient terms that describe the interface between the shear bands. The first is 
in the viscoelastic constitutive equation, with a characteristic (mesh) length I. The second is in the 
(generalised) Cahn-Hilliard equation, with the characteristic length ^ for equilibrium concentration- 
fluctuations. We show that the phase diagrams (and so also the flow curves) depend on the ratio 
r = l/^, with loss of unique state selection at r = 0. We also give results for the full shear-banded 
profiles, and study the divergence of the interfacial width (relative to I and ^) at the critical point. 

PACS numbers: 47.50.-|-d Non-Newtonian fluid flows- 47.20.-k Hydrodynamic stability— 36.20.-r Macro- 
molecules and polymer molecules 



I. INTRODUCTION 

For many complex fluids, the intrinsic constitutive 
curve of shear stress E as a function of shear rate 7 is 
non-monotonic, admitting multiple values of shear rate 
at common stress. For semi-dilute wormlike micelles, the- 
ory [1 H Q predicts the form ACEG of Fig. [T] In the 
range 7ci < 7 < 7c2 where the stress is decreasing, steady 
homogeneous flow (Fig. Ell) is unstable 01 . For an ap- 
plied shear rate 7 in this unstable range, Spenley, Gates 
and McLeish Q proposed that the system separates into 
high and low shear rate bands (711 and 7^; Fig. I^ls) and 
that any change in the applied shear rate then merely ad- 
justs the relative fraction of the bands, while the stress 
Ssci (which is common to both) remains constant. The 
steady state flow curve then has the form ABFG. Several 
constitutive models augmented with interfacial gradient 
terms have captured this behaviour 0, IE S 111 • 

Experimentally, this scenario has been widely ob- 
served in semidilute wormlike micelles inaini. The 

steady state flow curve (which is often attained only 
after very long transients has a well defined, re- 

producible stress plateau Egei. Goexistence of high and 
low viscosity bands has been observed by NMR spec- 
troscopy 10, 12, Hill- 

Further evidence comes from 
small angle neutron scattering (SANS) §|S|S|l3,|l3; 
and from flow birefringence (FB) |19l. i20l l2ll l2a|. which 
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FIG. I: Schematic flow curves for wormlike micelles: the 
homogeneous constitutive curve is ACEG; the steady shear- 
banded flow curve is BF (without concentration coupling in 
planar shear) or B'F' (with concentration coupling, or in a 
cylindrical Couette device). 



reveals a (quasi) nematic birefringence band coexisting 
with an isotropic one (but see (2.1 .24|). 

In some systems, the coexistence plateau is not per- 
fectly flat, but slopes upward slightly with increa sing 
shear rate (B'F' in Fig. See, for example, Ref. [23 
for GTAB(0.3M)/NaNO3(1.79M)/H2O at micellar vol- 
ume fraction (j) — 11%. This effect is much more pro- 
nounced in other, more concentrated systems that are 
near an underlying (zero-shear) isotropic-ncmatic (I-N) 
transition (0 « 30%) IHIllllg. 

In a cylindrical Gouette geometry, this upward slope 
is qualitatively consistent with the inhomogeneous stress 
arising from the cell curvature: as the high-shear band 
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FIG. 2: (a) Homogeneous shear rate and (b) banded profiles. 



at the inner cylinder expands outward with increasing 
apphed shear rate, the apphed torque must increase to 
ensure that the interface between the bands stays at the 
selected stress Egei ■ However a more generic explana- 
tion, independent of geometry, is that the shear banding 
transition is coupled to concentration 0,|23- In this case, 
the properties of each phase must change as the applied 
shear rate is tracked through the coexistence regime, be- 
cause material is redistributed between the bands as the 
high shear band grows to fill the gap. 

Generically, one expects flow to be coupled to concen- 
tration in viscoelastic solutions where the different con- 
stituents (polymer and solvent) have widely separated 
relaxation timescales jMSlMEillSlHllHllM- 
This was explained by Helfand and Fredrickson (HF) |23| 
as follows. In a sheared solution, the parts of an ex- 
tended polymer molecule (micelle for our purposes) in 
regions of lower viscosity will, upon relaxing to equilib- 
rium, move more than the parts mired in a region of 
high viscosity and concentration. A relaxing molecule 
therefore on average moves towards the higher concen- 
tration region. This provides a positive feedback mecha- 
nism whereby micelles can move up their own concentra- 
tion gradient, and leads to flow-enhanced concentration 
fluctuations perpendicular to the shear compression axis. 
This was observed in steadily sheared polymer solutions 
in the early 1990's [s^- In a remarkable paper, Schmitt et 
al. '2?! discussed the implications of this feedback mecha- 
nism for the onset of flow instabilities. Strongly enhanced 
concentration fluctuations were subsequently observed in 
the early time kinetics of the shear banding instability in 
Ref. [13 . 

Recently, therefore, we introduced a model of 
concentration-coupled shear banding |37l l38j by combin- 
ing the diffusive Johnson Segalman (d-JS) model 0, [s^ 
with a two-fluid approach (28L (iol l4ll li^ to concentra- 
tion fluctuations. This "d-JS-(/)" model does not address 
the microscopies of any particular viscoelastic system, 
but instead should be regarded as a minimal model that 
combines (i) a constitutive curve like that of scmidilutc 
wormlike micelles (Fig. ^ with (ii) the non-local (intcr- 
facial) terms required for selection of a unique banded 
state 6:] and (iii) a simple approach to concentration cou- 
pling- 

In Refs. |33,|33| I examined the linear stability of ini- 
tially homogeneous shear states in this d-JS-(/) model with 
respect to coupled fluctuations in shear rate 7, micellar 
strain W_ and concentration (f>. We thereby calculated 



the "spinodal" , inside which such homogeneous states are 
unstable. We also calculated the selected length scale at 
which inhomogeneity first emerges during startup flows 
in the unstable region. In the limit of zero concentra- 
tion coupling, the unstable region coincides with that of 
negative slope in the homogeneous constitutive curve, as 
expected; but no length scale is selected during startup. 
Concentration coupling enhances this instability at short 
length scales. It thereby broadens the region of insta- 
bility, and selects a length scale at which inhomogeneity 
must emerge. 

In the present paper, we compute the correspond- 
ing steady-state flow phase diagram (the "binodals" and 
their tie-lines). As far as we are aware, this is the 
flrst concrete calculation aimed at qualitatively describ- 
ing concentration-coupled shear banding for systems such 
as semi-dilute wormlike micelles. We start in Sec. ^ 
by describing the experimental background in more de- 
tail. We also compare our present calculation with the 
only other existing one for concentration-coupled shear 
banded states, in concentrated solutions of rigid rods [43j| . 
In Sec, mil we summarize our d-JS-(/) model. We then re- 
view our results for the spinodal onset of instability in 
Sec. IIVI In Sec. E] we describe our numerical procedure 
for computing the banded steady states, with brief dis- 
cussion of our careful study of mesh and finite size ef- 
fects. We then (Sec. IVI|I present our results for the fiow 
phase diagrams and shear-banded profiles. We conclude 
in Scclvni 



II. EXPERIMENTAL BACKGROUND; 
THEORETICAL CONTEXT 



In this section, we discuss in more detail the experi- 
mental evidence for concentration coupling in the shear 
banding of wormlike micelles. We survey both (i) concen- 
trated systems near the zero-shear I-N phase transition 
and (ii) semidilute systems, in which underlying nematic 
interactions are likely to be less important. Correspond- 
ingly, we compare the present calculation (aimed at the 
semi-dilute systems) with an earlier calculation of flow 
phase diagrams in concentrated solutions of rigid rods 
(near the I-N transition) j43] . 

The earliest observations of an upwardly sloping stress 
plateau in wormlike micelles were made by Schmitt et al. 

and Berret et al. 15, 26]. Schmitt et al. [IGj] studied 
CpClO3/NaClO3(0.05M)/H2O at the high miceflar vol- 
ume fraction (f) w 31%, just below the onset of the I-N 
transition at (t>i — 34%. In the steady-state flow curve, 
the stress increased smoothly up to the critical shear rate 
7^, where it showed a pronounced downward kink before 
curving upward again for 7 > 7£ (qualitatively like B'F' 
in Fig. SANS measurements confirmed a superposi- 
tion of nematic and isotropic contributions in this regime 
7 > 7f , with the nematic contribution rising linearly from 
zero at 7 = 7£. 

Berret et al. [H studied CpCl/hexanol/NaCl for 
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several micellar volume fractions, again at a volume frac- 
tion just below the onset of the zero-shear I-N transition 
« 32%). The overall height of the coexistence plateau 
(which again sloped upwards in 7) was found to fall with 
increasing surfactant concentration </> — > extrapolat- 
ing to zero at (/) ~ which is already biphasic in zero 
shear. They also found an increasing nematic contribu- 
tion to SANS patterns for increasing shear rates above 
7f. They further used the SANS data to show that the 
nematic (high shear) band was more concentrated than 
the low shear band. 

As noted above, the majority of existing calculations 
of shear-banded states have assumed uniform concentra- 
tion. An important exception is the calculation of Olm- 
sted and Lu Although this model was aimed at 
concentrated solutions of rigid rods, it broadly captured 
some of the experimental phenomenolo gy f or the concen- 
trated {(/) ~ 30%) wormlike micelles jlM IT^ |2^ . For 
example, the overall height of the coexistence plateau 
SsgI increased from zero as the concentration was reduced 
below the threshold (pi of the zero-shear I-N biphasic 
regime. The coexistence plateau sloped upward markedly 
in shear rate. In further agreement with experiment, the 
high shear (nematic) phase had a higher volume fraction 
of rods. It should be noted that model of Ref. was 
explicitly aimed at concentrated systems, which in zero 
shear are already close to the I-N transition: hence, the 
dynamics of the relevant order parameter Q was driven 

by a free energy that already contained a phase tran- 
sition. In contrast, the simple free energy FJW) we 
consider below has no underlying phase transition, and 
flow-induced efects are driven by convective, rather than 
dissipative (relaxational) dynamics. 
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FIG. 3: Height of the coexistence plateau in the system 
CpCI/NaSal/brine. The 5 leftmost points are taken from the 
data of Ref. 1^. The righthand point represents the zero- 
shear biphasic regime of this system, and is in accordance 
with the extrapolation of G{(j>) /TjbcI in Ref. (see main text 
for details). 

There have also been several experimental studies 



of concentration dependence in the shear banding of 
more dilute wormlike micellar solutions. For example, 
Berret et al. investigated the non-linear rheology of 
CpCl/NaSal/brine in the concentration range 5% — 20%, 
well below the I-Nc transition at « 36%. (Nc is "ne- 
matic calamitic".) In contrast to the more concentrated 
"prenematic" systems, the plateau height Sgoi decreased 
with decreasing concentration: see the left 5 points in 
Fig. 13 The width of the plateau also decreased so that 
the difference 7h — 7^ fell to zero at a critical point 
^c, 7c, Sc (leftmost p oint in Fig.l^J. (These trends are the 
same as those in Fig. 6(d) below.) In contrast the scaled 
plateau height Ssei/G'(0) (where G is the plateau modu- 
lus) decreased with increasing concentration and extrap- 
olated to zero in the zero-shear biphasic (I-N) regime at 
(j) ^ Ki 33%. According to this extrapolation (which 
is actually well beyond the final data point at « 22%), 
and in the absence of a divergence in G{(j)), the unsealed 
plateau height Sgoi must itself fall to zero at w 33% 
(rightmost data point in Fig. |3Jl , consistent with the be- 
haviour of the concentrated systems discussed above. To 
summarise, the plateau height Sgoi appears to be a non- 
monotonic function of concentration, increasing with (j) 
through the studied regime 5% < (j) < 20% before (prob- 
ably) falling to zero in the zero-shear biphasic I — N 
regime (f> « 33%. 

Although this experiment showed that shear band- 
ing depends on the overall concentration of the solu- 
tion, there is relatively little evidence for concentration- 
coupling [i.e. concentration differences between the 
bands) in such dilute systems, far from the I-N tran- 
sition. Indeed, the experiment just described revealed 
no discernible upward slope in the coexistence plateau. 
We are not aware of any measurements of concen- 
tration differences between the coexisting shear bands 
in such systems. Nonetheless, recent experiments on 
CTAB(0.3M)/NaNO3(1.79M)/H2O at </) = 11% jH did 
reveal a stress plateau with slight upward slope. Along 
with the generic expectation that flow should be cou- 
pled to concentration in these viscoelastic solutions, this 
suggests that an explicit calculation of concentration dif- 
ference in the shear bands of systems far from an I-N 
transition might be worthwhile. 

In this paper, therefore, we present the first such cal- 
culation, using our d-JS-0 model [s^, HI]- In contrast 
to the work of Olmsted et al. ^3 fo^' concentrated rigid 
rods, shear banding in the d-JS-(/) model is not due to 
any underlying nematic feature of the elastic free energy 
F'^(W_). Instead the instability results mainly from the 
non-linear effects of shear (the intrinsic constitutive curve 
has a region of negative slope), though it can be strongly 
enhanced by concentration coupling in systems close to 
an underlying Cahn-Hilliard (CH) demixing instability 
(governed by the osmotic free energy F° (</))). Indeed, 
the d-JS-(/) model captures a broad crossover between 
(i) instabilities that are mainly mechanical (governed by 
the negative slope of the flow curve) and (ii) instabilities 
that are essentially CH demixing (governed by F° {(/))), 
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but now triggered by shear. [Likewise, in practice there 
should be no sharp distinction between concentrated mi- 
celles with an underlying nematic feature in F'^{W_) on 
the one hand and "non-nematic" (more dilute) systems 
on the other: any more refined model should allow a 
smooth crossover between the two_ cases. This will be 
the focus of a future publication (44i|.] 



III. MODEL 

In this section we outline the d- JS-(/) model, which cou- 
ples shear banding instabilities to concentration in a sim- 
ple way by combining the non-local Johnson-Segalman 
(d-JS) model 39] with a 2-fluid framework (H l3l| for 
concentration fluctuations. While this description is self- 
contained, readers are referred to Ref. for fuller de- 
tails. 



A. Free energy 

In a sheared fluid, one cannot strictly define a free en- 
ergy. Nonetheless, for realistic shear rates, many internal 
degrees of freedom of a polymeric solution relax quickly 
on the timescale of the moving constraints and are there- 
fore essentially equilibrated. Integrating over these fast 
variables, one obtains a free energy for a given fixed con- 
figuration of the slow variables. For our purposes, the rel- 
evant slow variables are the fluid momentum and micellar 
concentration (j> (which are both conserved and therefore 
truly slow in the hydrodynamic sense), and the micel- 
lar strain W_ that would have to be reversed in order to 
relax the micellar stress (which is slow for all practical 



purposes): 



dR' dR'p 



(3.1) 



where SR' is the deformed vector corresponding to the 
undeformed vector SR. 

The resulting free energy is assumed to comprise sep- 
arate osmotic and elastic components, 



F = F°((/)) -f i^'=(^. 
The osmotic component is 

F°(0) = 



(3.2) 



/(0) + |(W 



d^(7(i + ev)/"i0(9)i 



(3.3) 



where /" is the osmotic susceptibility and ^ is the equi- 
librium correlation length for concentration fluctuations. 
The elastic component is 

F^'i.K, = 5 / [K - ^°g(| + W] (3-4) 

in which G(0) is the micellar stretching modulus. 
B. Dynamics 



The basic assumption of the two-fluid model is a sep- 
arate force balance for the micelles (velocity v^^; volume 
fraction (f>) and the solvent (velocity v^) within any ele- 
ment of solution. 

The micellar force balance equation is: 



Pm idt+v^.V)v^ = Y.G{^)W^ ^^^^ + 2 V. (/.r7m^° - C(0)l^.ci - (3-5) 

Q(p m 

In this equation, G{(j))W_ = '2'G{4>)W_.j^ is the viscoelastic micellar backbone stress due to deformation of the local 
molecular strain, while the osmotic stress results from direct monomeric interaction. The Newtonian stress 
2077nii2'^^ describes fast micellar processes {e.g. Rouse modes) with if'^^ the traceless symmetric micellar strain 
rate tensor. The force Cl^rd, where w — y_^^ ~ v^, impedes relative motion; ^ is the drag coefficient (Eq. I3.8|l . 
Incompressibility determines the pressure p. 

Likewise, the solvent force balance comprises the Newtonian viscous stress, the drag force (equal and opposite to 
the drag on the micelles), and the hydrostatic pressure: 

p,(l - 0) {dt + v,.V)v, = 2V. (1 - cb)Ti,Rl + C(0)2;,ei - (1 - 0)Vp. (3.6) 

Equations 13 . 51 and 13 . 61 contain the basic assumption of "dynamical asymmetry", i.e. that the viscoelastic stress acts 
only on the micelles and not on the solvent. Adding them, and assuming equal mass densities pm ~ ps ^ p, we obtain 
the overall force balance equation for the centre of mass velocity, v = (f>v^ + {l — (/))Ws: 

p(5t-hz).V)w = Aw=Y.G(0)li:-0V^^^-l-2V.(/)7y,„Z2° + 2V. (1 - 0) ^ Vp. (3.7) 
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Subtracting them (with each predivided by its own volume fraction), we obtain an expression for the relative motion 



- Vg, which in turn specifies the concentration fluctuations: 



= -Y • 0(1 - (b)v,,i = -V • 



0^(1-0)^ 
C(0) 



V-G{(I>)K ^SF (j)r]mRl^ 2V • (1 - </>) r/, £° 



(3.8) 



which defines the micellar diffusion coefficient D = f"{(i))<t>^ (1 — 0)^/C- We have omitted negligible inertial corrections 
to Eqs. 1371) and W^. 

The essence of the 2-fluid model is that the physically distinct elastic and osmotic stresses appear together in the 
force-balance equation H3.7|l and also in the generalised CH equation (|3.8() . This allows micellar diffusion in response 
to gradients in concentration and in the viscoelastic stress. We will see below that this gives rise to a positive HF 
feedback between concentration and flow [S^, allowing micelles to diffuse up their own concentration gradient. 

For the dynamics of the viscoelastic micellar backbone strain we use the phenomenological d-JS model [5l.l39|: 



{dt + v^-V)W_ = a[D .W + W.D ) + {W.n .W) + 2D 



W 



r(0) Ti(b) 



(3.9) 



where 2^^^ = Yw^-(VWjjj)^ with (Vw^Ja/j = da{vm)f3- 
t{4)) is the Maxwell time and I is a length scale discussed 
in Sec. IIII Dl below. The slip parameter a measures the 
non-affinity of the molecular deformation, i.e. the frac- 
tional stretch of the polymeric material with respect to 
that of the flow field. For \a\ < 1 (slip) the intrinsic 
constitutive curve in planar shear is capable of the non- 
monotonicity of Fig. 

We use Eqns. EI| EH and EiH together with the in- 
compressibility condition, V_.v = 0, as our model for the 
remainder of the paper. 

C. Flow geometry. Boundary conditions 

We consider idealised planar shear bounded by infinite 
plates at y = {0,L} with (w,Vu,V Au) in the {x,y,z) 
directions. We allow variations only in the flow-gradient 
direction, and therefore set all other derivatives to zero: 
dx . ■ ■ = 0, dz . . . = 0. In appendix 1X1 we give all the rele- 
vant components of the model equations 13. 71 13.81 and 13.91 
in this coordinate system. 

The boundary conditions at the plates are as follows. 
For the velocity we assume there is no slip. For the con- 
centration we assume 

dycj) = dl(j) = 0, (3.10) 

which ensures (in zero shear at least) zero flux of con- 
centration at the boundaries. Following Ref. |5(|, for the 
micellar strain we assume 

dyWaf3 = 0y a,(3. (3.11) 

Conditions 13 . 101 and 13 . 1 II together ensure zero concentra- 
tion flux at the boundary even in shear. For the con- 
trolled shear rate conditions assumed throughout, 

- /"^ 

j — dyj{y) — constant. (3.12) 

^0 



D. The interfacial terms 



The model contains two different interfacial terms. 
The first is the gradient term on the RHS of Eqn. (|3.9|l . 
The length I in this term could, for example, be set by 
the mesh size or by the equilibrium correlation length for 
concentration fiuctuations. Here we assume the former, 
since the dynamics of the micellar conformation are more 
likely to depend on gradients in molecular conformation 
than in concentration. Physically, one can interpret the 
gradient term in equation 13.91 as resulting dynamically, 
from the diffusion of stretched molecules across the inter- 
face [45I I , or statically, from nematic interactions between 
the micelles, or both. There is, at present, no accepted 
theory for these gradient terms in semi-dilute solutions. 
The equilibrium correlation length ^ of course still enters 
our analysis through our second interfacial term, in the 
osmotic free energy of Eqn. 13.31 

Together, I and ^ set the length scale of any interfaces. 
Throughout this paper, we study the physical limit in 
which / and ^ are small compared to the system size so 
that we have a sharp interface connecting two bulk homo- 
geneous phases. In this case, the solution to Fans. 1X71 13. 81 
and 13.91 naturally fits the zero-gradient boundary con- 
ditions, and is invariant under y — > y/2, I 1/2 and 
^ — > f/2. Therefore, a simultaneous reduction in I and ^ 
by the same factor only changes the overall length of the 
interface, and not the values of the order parameters in 
each phase (which determines the phase diagram). How- 
ever the phase diagram does depend slightly on the ratio 
r = l/^: below we will give results for r = 0, r = 00 and 
r — 0(1). This provides a concrete example of the early 
insight of Lu and co-workers that the banded state 
must depend on the nature of the interfacial terms. This 
contrasts notably with equilibrium phase coexistence, in 
which the dynamical equations are integrable and there- 
fore insensitive to interfaces. 
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Parameter 


Symbol Q 


Value at = 0.11 


dlogQ 
d log (f) 


Rheometer gap 


L 


0.15 mm 





Maxwell time 


r 


0.17s 


1.1 


Plateau modulus 


G 


232 Pa 


2.2 


Density 


P 


lO^kgm"^ 





Solvent viscosity 




lO^^kgm^^s"^ 





Rouse viscosity 


Vm 


0.4kgm"^s'^ 





Mesh size 


I 


2.6xlO"V 


-0.73 


Diffusion coefficient 


D 


3.5 X lO-^m^s"^ 


0.77 


Drag coefficient 


c 


2.4 X lO^^kgm-^s-i 


1.54 


Correlation length 




6.0 X 10"^m 


-0.77 


Slip parameter 


a 


0.92 






TABLE I: Experimental values of the model's parameters at 
volume fraction (p = 0.11 (column 3). Scaling laws for the de- 
pendence of each parameter upon (j> (column 4). In most cal- 
culations we use the reference values of column 3 at = 0.11, 
then tune <j) using the scaling laws of column 4. Only where 
stated do we allow the parameters to vary independently. 



E. Model parameters 

The d-JS-0 model (Eqns. EZ| EHI EHll has the follow- 
ing parameters: the solvent viscosity rjs and density p; 
the plateau modulus G; the Maxwell time t; the Rouse 
viscosity 77 m; the mesh size I; the osmotic modulus /"(0) 
and the equilibrium correlation length ^ (recall Ean. l3.3ll : 
the drag coefhcient C and the slip parameter a. We also 
need to know the typical rheometer gap, L. A reference 
set of parameter values at </> = 0.11 is summarised in 
table m These values were taken from experiment or cal- 
culated using scaling arguments: see Ref. [s^ for details. 
Note that explicit data is not available for f"{4>); however 
dynamic light scattering gives the diffusion coefficient 



D = 



r 



^(1 



(3.13) 



In this paper we will be guided by these parameter values, 
but subject to the following considerations. 

First, we are only interested in steady states so for 
convenience can take the limit of zero Reynolds number 
(p = 0) and rescale the kinetic coefficient 1/C so that 
the diffusion time L^/D is of order the Maxwell time. 
These choices have no effect on the steady state, but 
make our numerical calculation of it much more efficient 
(by evolving the dynamical equations 13.71 13.81 and I3.9|l . 
Second, realistic interfaces are much narrower than the 
typical rheometer gap, with I and 1^ both of O{10~*L). 
To resolve such interfaces (allowing a minimal 10 numer- 
ical mesh points per interface) would therefore require 
O(IO^) grid points, while in practice we are limited to 
O(IO^). We will therefore use artificially large values of 
I and ^. However this does not affect the phase diagram, 
provided the interface is still small compared with the 




FIG. 4: Intrinsic flow curves (dotted lines) for (j> = 
0.11, 0.10 ... 0.01 (downwards). Spinodals for concentration 
couplings a = 10"^, 10~^ 10"". 



gap size: see Sec. El for details, and Fig. [Slfor a typical 
banded profile. Finally, we artificially increase the Rouse 
viscosity 77 m by a factor 50 to ensure, again for numerical 
convenience, that the shear rate of the high shear phase is 
not too large. This does quantitatively change the phase 
diagram, but we checked that the qualitative trends are 
not affected. 

Exploring this large parameter space is a daunting 
prospect so we shall not, in general, vary the parame- 
ters independently of each other. Instead we simply tune 
the concentration 0, relying on known semi-dilute scaling 
laws for the 0-dependence of the other parameters (col- 
umn 4 of tableP) . However we will, in separate sweeps, 
vary the degree of concentration coupling, which is dic- 
tated by ratio of the elastic term ^■G{(j))W_ to the osmotic 
term and which we encode in the parameter 



G"(0 = O.ll) 



2/"( 



0.11) 



(3.14) 



(where a prime denotes a derivative) . In other ^—sweeps 
we will vary the characteristic interface widths Z(0 = 
0.11) and ^(0 = 0.11), to investigate any dependence 
of the phase diagram on the ratio l/^ in the double limit 
1/ L ^ 0, L ^ 0. In what follows, we adopt the conve- 
nient shorthand of Z for ?(</) = 0.11) with the understand- 
ing that I does actually varies with (j) according to the 
scaling given in table HJ We do likewise for ^. 

Throughout we rescale stress, time and length so that 
G{(j) = 0.11) = 1, T{(j) = 0.11) = 1, and L = 1. 



IV. INTRINSIC FLOW CURVES; SPINODALS 

The homogeneous intrinsic steady state flow curves 
i;(7,0) = G{4>)Wxy + 77(^)7 that satisfy dtv = dt4> = 
dtW = are shown as dotted lines in Fig. 21 (The av- 
erage viscosity ri{(j)) = (t^Vm + (1 ~ 4')ls-) The region 
of negative slope ends at a "critical" point 0c ~ 0.015. 
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CPCl/NaSal in brine Q shows the same trend. For com- 
pleteness, in App. we give analytical results for the 
steady state conditions in homogeneous shear flow. 

InRef. wc linearised in fluctuations about these 

homogeneous states to find the spinodal region in which 
the homogeneous states are unstable. The spinodals are 
shown in Fig.01for different levels concentration coupling, 
a. 

In the limit of zero concentration coupling a — > 0, 
fluctuations in the "mechanical variables" , W_ and 7 de- 
couple from those in concentration, and are unstable in 
the region of negative constitutive slope, as expected. 
Separately, the concentration could have its own Cahn- 
Hilliard demixing instability, when the diffusion coeffi- 
cient D < 0; however we are interested only in fiow- 
induced instabilities and set D > throughout. For finite 
a > 0, the region of mechanical instability is broadened 
by coupling to the concentration fluctuations, as seen in 
Fig. ^ This can be understood as follows. Consider the 
first term in the square brackets of Eqn. (|3.8(l . This causes 
micelles to move up gradients in the viscoelastic stress W_, 
thereby increasing the concentration in stressed regions. 
If G'{4>) > (assumed here), the increased concentration 
causes the stress to increase further, closing a positive 
HF feedback loop whereby the micelles can diffuse 
up their own concentration gradient. 



V. NUMERICAL DETAILS 

In this section, we outline our numerical procedure for 
solving the dynamical eg nations 13.71 1X^1 and 13 . 9l and dis- 
cuss our careful study of time-step, mesh size and finite 
size effects. Readers who are not interested in these is- 
sues can skip this section. 

We consider variations only in the flow gradient direc- 
tion, in which we discretise y e 0, 1 on an algebraic grid 
Vn — n/Ny for n = 0, 1 . . . Ny. We stored (j) and W_ on 
these grid points. The velocities and were stored 
on half grid points ?/„+i/2, and we used linear interpola- 
tion between the half and full grid points. Likewise we 
discretized time such that t„ — n/S.t. We evolved the 
discretized equations 13.71 13.81 and 13.91 using the Crank- 
Nicholson algorithm which is semi-implicit in time, with 
centred space derivatives 

For each run, we seeded an initial profile that was ei- 
ther homogeneous up to a small random contribution, 
or inhomogeneous according to (p — + Acos(7ra;)] 
(with A K, 0.1). We then evolved the discretized equa- 
tions 13 . 71 13 . 81 and l3 . 9l under an imposed wall velocity until 
a steady banded state was reached. We checked that the 
homogeneous phases between the interfaces were insen- 
sitive to the initial conditions. However, for the random 
initial condition several bands could form (and did not 
coarsen over any accessible timescale) . Therefore in most 
runs we used the co-sinusoidal initial profile, to conve- 
niently obtain just two bands (as in Fig. 112(1 . 

For the dynamics to be independent of time-step. 
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FIG. 5: Main figure: steady banded concentration profile at 
a = 10"^ ^ = 0, i- = 7.0, 4> = 0.15. Solid lines / = 0.016, 
for iNy,At) = (100, 0.05), (200, 0.0125), (400, 0.0125), dotted 
lines I = 0.008 for iNy,At) = (200, 0.05), (400, 0.003125), 
dashed lines I = 0.004 for {Ny,At) = (400, 0.05), (800, 0.05). 
Upper inset, the same data, enlarged in the left hand phase 
(decreasing with increasing Ny). Lower inset: correspond- 
ing selected stresses, for the same parameter values and mesh 
sizes (decreasing stress with increasing Ny). 



a very small time-step has to be used. However the 
numerically-attained steady is much less sensitive hence 
allowing much larger time-steps. A typical steady state 
changes by less than 10^'^ % for a factor-two reduction 
in timestep El For the special case of ^ = 0, timesteps 
At oc Ny'^ can be used, since the highest spatial deriva- 
tive is second order. For ^ 7^ 0, we have a fourth or- 
der derivative in Eqn. 13.81 and much smaller timesteps 
At cx Ny'^ must be used. 

In all our calculations, we are interested in the physi- 
cal limit where the interface width is much smaller than 
the rheometer gap. This creates a delicate balance, since 
narrow interfaces require a very fine grid. Therefore we 
adopted the following procedure. For any fixed value of 
the interfacial lengthscales I and ^, we performed several 
runs with progressively finer meshes (but always with a 
small enough time-step) until the shear banded profile 
and selected stress didn't depend on the mesh. This is 
quite easy to achieve: a typical steady state presented be- 
low changes by less than 0.1 % upon doubling the number 
of mesh points. We then reduced I and ^ (in fixed ratio) 
until the order parameters in the homogeneous phases 
changed by less than 0.5% upon further halving of / and 
^ (but always ensuring convergence with respect to the 
number of grid points) . A sample study of these issues is 
presented in Fig.|5lfor the special case ^ = 0. 



VI. RESULTS 

We now present our results for the steady-state flow 
phase diagrams, flow curves and shear banded profiles. 
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(a) Tie lines (solid); spinodal (dashed). 



(b) Tie lines (solid); spinodal (dashed); homogeneous 
constitutive curve (dotted); for I = 0.008. 
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(c) Tie lines (solid); spinodal (dashed) for I = 0.008. 



(d) Coexistence regime of the macroscopic flow curves 
(solid); spinodal (dashed); homogeneous constitutive 
curves (dotted) 



FIG. 6: Phase diagrams and flow curves for a = 10~^, ^ = 0.0 for small l/L. (Recall that / is actually a function of (f): we are using 
the convenient shorthand of / for the value [{(f} — 0.11).) (a) Thin (upper) solid lines: tie lines for / = 0.016, Ny — 100, At = 0.05. 
Thick (lower) solid lines: tie lines for / = 0.008, Ny = 200, At — 0.05. As described in the main text, we actually rescaled I 
in the successive runs of each 0— sweep (i.e. as <j) was tracked from 0.15 down to <^c) so that the interfacial width remained 
(approximately) constant throughout the sweep: the value of I in the figure legends refers to the value used in the first run of the 
sweep, at = 0.15. (b,c) Solid lines: tie lines repeated in the (S,7), (S,</)) representations for I = 0.008, A^j, = 200, At = 0.05. 
(d) Solid lines: macroscopic flow curves for (j) = 0.11, 0.10, . . . 0.04 (downward). These flow curves were recontructed from the tie 
lines of the phase diagrams (using the tie lines shown in this figure, and some additional ones). Because we have only calculated 
tie lines for discrete values of E, in some cases the reconstructed fiow curves stop short of the single-phase region, and have 
been continued by eye with a dashed line. The inset in (d) shows the same data, but on a log-log plot. The experimentally 
observed slope 0.3 is marked as a dot-dashed line for comparison. The spinodal is shown in each of Figs a-d as a dashed line. 
In (b,d) the thin dotted lines are the intrinsic (homogeneous) constitutive curves for (f> = 0.11, 0.1 .. . 0.01 (downwards). 



Because one of our aims is to show that the shear- 
banded state depends on the nature of the interfacial 
terms, we consider three separate cases: (A) interfacial 
terms only in the viscoelastic constitutive equation 13.91 



{I ^ 0, £^ ^ 0, r = I /(, — oo); (B) interfacial terms in both 
the constitutive and concentration equations . 13 . 91 and 13.81 
/ 0,^ 7^ 0, r = 0(1)); and (C) interfacial terms only 
in the concentration eauation l3.8l = 0, ^ 7^ 0, r = 0). 
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A. Interfacial terms only in the viscoelastic 
constitutive equation: 1^0, ^ = 0. 

In this section, we set the correlation length for con- 
centration fluctuations, ^, to zero and consider small but 
non-zero values the interfacial lengthscale I in the consti- 
tutive equation 1^21 



1. Flow phase diagrams 
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FIG. 7: Phase diagrams for three different degrees of coupling 
to concentration for ^ = and small l/L. 

For any given shear banded profile, the values of the 
order parameters in each of the two homogeneous phases 
specify the two ends of one tie line in the phase diagram. 
Analogously to equilibrium tie lines, the concentrations 
and strain rates of the coexisting states are related to 
the mean strain rate 7 and mean concentration (jj by the 
lever rule, 



^ = /3(/.i + (1 - (3)^2 
^ = /37i + (1-/3)72, 



(6.1) 
(6.2) 



where P is the volume fraction of material in state 
{4>i , 7i ) • For each of several values of the concentra- 
tion coupling, a, we calculated the full phase diagram 
via a succession of shear startup runs, all at the critical 
shear rate 7c (a) (determined from Fig. for average 
concentrations ranging from </) = 0.15 down to the crit- 
ical value 0c (a)- For concentrations below the critical 
point the response of the system is smooth as a func- 
tion of stress. In our model, this arises because decreas- 
ing concentration reduces the viscosity of the low shear 
rate branch faster than it reduces the viscosity of high 
shear rate branch. Hence the stress maximum decreases 
with decreasing concentration, disappearing when the 
stress maximum vanishes. Alternatively, in a more di- 
lute system the plateau modulus and Maxwell time are 
both smaller, and one expects a smaller stress and higher 
strain rate at the onset of instability. 

The results for a = 10^^, which gives rather strong 
concentration coupling, are shown in Fig.Et,b,c. Because 
the width, S, of the interface in the banded state is set by 



Z, but with a prefactor that diverges at the critical point, 
in each successive run we rescaled I so that S remained 
(approximately) equal to its value (^ L) in the first run 
at (f> ~ 0.15. We return below to study the divergence of 
6/1 at the critical point. 

To illustrate the finite size considerations of Sec. |V] 
(above), in Fig. 6(a) we show the tie lines obtained for 
two different (starting) values of I. All the results are 
converged with respect to mesh fineness and timestep 
(not explicitly shown), but the tie lines differ slightly 
between the two values of /. However all seem to be 
consistent with one given binodal line: we do not have 
any explanation for this apparent consistency. 

To investigate the effect of reducing the coupling to 
concentration, we repeat the phase diagram for a — 10^^ 
alongside that for a = 1Q~^ and a = lO"** in Fig.[7| As 
expected, the concentration difference between the bands 
tends to zero as a — > 0. 



2. Flow curves 

So far, we have discussed the flow phase diagrams. 
Measurement of these diagrams still presents an open 
challenge to experimentalists, due to the difficulty in 
measuring the concentration of micelles in each band (al- 
though SANS data has been used to estimate the bands' 
concentrations in systems near the I-N transition 26]). 
In this section we discuss the macroscopic flow curves, 
which are relativly easily measured using conventional 
bulk rheology. However it is important to realise that a 
s_et of flow curves I](7,0) measured for several values of 
(j) actually contains the same information as the phase 
diagram: reconstruction of the latter from the former is 
described in Fig.^|A full set of flow curves could there- 
fore be used to check measurements of concentration dif- 
ferences. 

In this work, we take the opposite approach for con- 
venience, and reconstruct the steady-state flow curves 
from the tie lines of the phase diagram. The results are 
shown in Fig. 6(d) The inset shows the same data on a 



log-log plot, to enable comparison with Ref. 15] in which 
the coexistence plateau in a log-log representation was a 
reasonably straight line (over the shear-rate range inves- 
tiga ted) w ith slope 0.3. Note that the results shown in 
Fig. |6(d)[ar e in units of G{(t> = 0.11) and r(</) = 0.11). 
In Ref. |46| . Berret replotted the flow curves in units 
of G{(j)) and r(0), finding scaling collapse of the family 
I](7, (/))/G(0) vs. 7t((/)) in the low shear regime 7 — > 0. 
We do not find this scaling collapse (Fig. Ill|) because we 
have used an artificially large high-shear Newtonian con- 
tribution 777 for numerical convenience (recall Sec. lIIIE|l : 
the overall zero shear viscosity, G(0)t(0) -1-77(0) therefore 
does not scale as G'((/))t((/)), even approximately. 

To check the reconstruction of flow curves from the 
phase diagram, we also explicitly calculated the flow 
curve at a single = 0.11. To do this, we first per- 
formed a shear startup at a given 7 in the unstable region. 
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FIG. 8: Macroscopic flow curves for a = 10"^ at = 0.11. 
Thick solid line: reconstructed from the tie lines of the phase 
diagram. Dot dashed and dashed lines: calculated by directly 
measuring the average stress and strain rate during a strain 
rate sweep for I = 0.016, Ny = 100, At = 0.05 (dot-dashed) 
and I = 0.008, Ny = 200, At = 0.05 (dashed). (The slight 
discrepancy between these three curves is discussed in the 
text.) The thin dotted line is the intrinsic (homogeneous) 
constitutive curve. 
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wards in 7. In Fig.|2|we compare the (directly measured) 
macroscopic flow curve for the three levels of concentra- 
tion coupling shown in Fig. [7| as expected, the slope of 
the flow curve tends to zero with the degree of concen- 
tration coupling a. 

The upturn in the measured flow curve at the edge of 
the coexistence plateau (apparent at the lower binodal 
for a = 10^^ in Fig.^ results again from the finite value 
of S/L: the interface bumps into the edge of the rheome- 
ter when one of the bands gets very narrow. We expect 
this (steady-state) effect to be much less pronounced in 
experimental systems, since realistic interfaces are much 
smaller than those used in our numerical study. Only 
near a critical point, where the interface becomes very 
broad (for fixed I), would we expect to see a true steady- 
state bump at the edge of the plateau. Nonetheless, pro- 
nounced bumps are often apparent in data obtained via 
upward strain-rate sweeps. However in most cases this is 
likely to be a metastable effect, so that the bump could 
be eliminated (or at least reduced) by reducing the rate 
of the sweep [Tl| . 

As noted in Sec. U in a curved Couette geometry the 
"plateau" (B'F' of Fig.P) in the flow curve will slope up- 
wards due to the inhomogeneity of the stress field, even 
without concentration coupling. It should be noted that 
all calculations in this paper are for a planar shear ge- 
ometry, and the slope of our flow curves in the coexis- 
tence regime results solely from concentration coupling. 



In fact, the slope in Fig. 



is far greater than one 



would typically expect from curvature effects: for a Cou- 
ette cell with radius R and gap 6R, the stress measured at 
the inner Couette wall would change by = 25R/R 

over the coexistence regime, and so too would the relative 
change in torque through the coexistence "plateau" . The 
slope of Fig. 6(d) would therefore require an atypically 
large curvature oi SR/R ^ 0.5. 



FIG. 9: Macroscopic flow curves (from direct measurements 

of the stress and strain rate) for three different degrees of 3. Interfacial profiles; divergence of interface width at the 
coupling to concentration. critical point 



We then (without reinitiahsing the system) decreased 7 
in steps to the edge of the coexistence regime, ensuring 
that a steady state was reached before measuring the to- 
tal stress. We then reinitialised the system and repeated 
the entire procedure, but now with increasing 7— jumps. 
The results are shown in |S1 for two different values of 
I. The slight discrepancy between the directly measured 
flow curve "plateaus" (i.e. the inhomogeneouse part of 
the flow curve) and those reconstructed from the tie lines 
is due to the finite size of the interface S relative to the 
cell L, and so is smaller for the smaller value of d/L. The 
construction described in Fig. llOl implicitlv assumes that 
S/L = 0. 

As expected for this value of a (which gives a large 
concentration difference between the bands; Fig. |6(a)1 l, 
the steady state flow curve "plateau" slopes strongly up- 



We now turn to the interfacial profiles and widths. A 
full steady state banded profile for a — 10"^ (correspond- 
ing to the rightmost/uppermost tie line in fig. |SJi,b,c) is 
shown by the thick lines in Fig. 1121 As required, the in- 
terface is smooth on the scale of the mesh, but sharp on 
the scale of the gap size, i.e. L/Ny <C (5 <C i = 1 where 5 
is the width of the interface. Note that the shear rate is 
negative across the gap since we have chosen to move the 
wall at ?/ = 0; accordingly we have plotted —Wxy, since 
Wxy is antisymmetric in shear rate. —W^y is rather small 
in the high shear band, as expected from the underlying 
constitutive non-monotonicity. Meanwhile Wxx is very 
large, while Wyy « —0.5 (recall that W_ measures defor- 
mation relative to the unit tensor 5): this corresponds 
to the micelles being highly stretched along the flow di- 
rection and is consistent with the experimental observa- 
tion that the first normal stress difference progressively 
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FIG. 10: Reconstruction of the flow phase diagram from a family of ma croscopic flow curves E(7, </!)), measured for several 
different average concentrations 4>. Consider the flow curves of Fig. [T0(a)1 The curve that starts at A and ends at B is for an 
average concentration (j) = 0.08. Points A and B are at the edge of the two-phase region. Reading off the stress from Fig. |10(a)1 



A a nd B g ive use two points on the binodal in Fig . |10(b ) 
Repeating this for all the circles in Fig. |10(a) 



Likewise reading off the strain rate, we get points A and B in 
we can construct many points on the binodal in Figs 



Fig. mHc) _ _^ 

|10(c)| which can then be interpolated over to give the full binodal. We now just need to specify the tie lines. In Fig, 



io(b: 



Wh 



and 
this 



is trivial: all tie lines are horizontal since the coexistence occurs at common stress (for gradient banding). In Fig. |10(c)| to get 
the slope of the tie line that star ts at B we proceed by recalling that the tie line represents constant shear stress. Therefore 
we find another point, D, in Fig. |10(a)1 that is at the same stress as point A, and read off its average strain-rate. Its average 
concentration is already known. This gives point C in Fig. |10(c)| Similarly, D is the image of point B at constant stress. 
Repeating this process we can fill in all the tie lines of the phase diagram. 



E/G((|)) 




10 
Yt:((^) 



FIG. 11: Macroscopic flow curves as shown in Fig. |6(d)] above, 
but now with the stress in units of G{(j)) and the strain rate 
in units of r (</>). 



increases throughout the banding regime [ig. The con- 
centration is lower in the high shear band, where Wyy 
is smaller (more negative): this is a direct result of the 
tendency of mielles to move up gradients in Wyy, as de- 
termined by Egn. 13. ijl above. 

In fact the interface width, (5, is slightly different for 
each order parameter: we define it to be the distance 
between the two points where the change in that order 
parameter between the two homogeneous phases is 25% 
and 75% complete. For a fixed value of I (which sets 
the overall scale of the interface width) , 5 diverges at the 
critical point (for each order parameter). In tracking 
down towards the critical point, therefore, we continually 
rescaled I to ensure that the interface width remained 



approximately constant. In each case, we measured 5/1, 
for each of Wxy, Wxx, Wyy and (j): see Fig.^l According 
to mean field theory, the divergence should be of the 
form 6 /I ^ (S — Sc)^^/^. The power —1/2 is accordingly 
shown in Fig. ^| as a guide for the eye. 



B. Interfacial terms in both the viscoelastic 
constitutive equation, and in the concentration 
equation: I ^ 0, ^ 0. 

We now study the effect of including interfacial gra- 
dient terms in the concentration equation 13.81 (so that 
now ^ ^ 0) as well as in the viscoelastic equation 13.91 
1^0. Hence, while in the previous section we consid- 
ered r = Ij^ = oo, then, we now consider r = 0(1). 
In Fig, 
Cor 
Fig 



14(a) 



we give the phase diagram for r = 0.4. 
Comparing it with our results for r — oo (also shown in 
I, we see that the slopes of the tie lines and the 
overall binodal both depend quantitatively on r. [The 
difference between the results for r = oo and r = 0.4 
is far greater than any "error" associated with the fact 
that we are not quite in the limit At — s- 0, INy — > oo, 
S^Ny ^ OO, //L — > and ^/L 0.] This provides a 
concrete example of the fact that shear-banding coexis- 
tence is determined by, and non-universal with respect 
to, the interfacial terms 6j. As noted above, this con- 
trasts sharply with the equilibrium case, in which the 
equations of motion are integrable and so the phase di- 
agram is independent of the interfacial terms. Although 
conceptually important, this dependence is in practice 
rather weak: the overall features of the phase diagram 
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FIG. 12: Steady state shear banded profile at q = 10 , 
7 = 4.64, (j> = 0.15 for two different ratios r = i/^. The 
thick lines are for r = oo = 0.008, C = 0.0), Ny = 200, 
= 0.05, as considered in this section. The thin lines show 
the corresponding resuHs for r = 0.4 (/ = 0.00 8, j = 0.002), 
Ny = 200, = 0.00625 (to be discussed in Sec. lVrSl below). 
for comparison. 
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FIG. 13: Scaled interface width 5/1 versus the distance from 
the critical stress S — Ec. The dotted line is a power —0.5. 



are unchanged. The critical point is unaffected. 



In Fig. 14(b) we show the corresponding macroscopic 



flow curves, reconstructed using the tie lines of Fig. 14(a)_ 
Because we only calculated a few tie lines in this case, 
recontruction is rather sparse. Nonetheless, the slight 
difference between r = oo and r — 0.4 is apparent. 

In Fig. El we compare a full banded profiles for r — 
oo and r = 0.4. The slight dependence on r is again 
apparent. 
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(b) Partially reconstructed macroscopic flow curves 
(solid); homogeneous constitutive curve (dotted). 



FIG. 14: Phase diagrams and flow curves for a — 10~^ and 
r = 0.4 (l = 0.008, ^ = 0.02) with the corresponding data for 
a = 10~^ and r = oo (l = 0.008, C = 0) for comparison. 



C. Interfacial gradient terms only in the 
concentration equation: / = 0, ^ 7^ 0. 

Finally we set the interfacial length I in the constitu- 
tive equation equal to zero. The constitutive equation is 
now local, and the only source of spatial gradients is the 
equilibrium correlation length for concentration ffuctua- 
tions (Bans. ITS! and in^ : r = l/^ = 0. In the absence of 
concentration coupling, a = 0, it is known that there is 
no uniquely selected, smoothly shear banded state when 
Z = I3 • Here we investigate whether a smoothly banded 
state is selected for a ^ 0, by virtue of the interfacial 
terms in the concentration equation. 

Our numerics only gave a smoothly banded profile for 
stresses near the critical point, even for the largest ac- 
cessible values of ^ and Ny . The profiles shown from left 
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(c) Steady-state profiles for initial conditions 
4>{y) = (j> + OAcos{iTy) and (f>{y) = if> + 0.7 cos{iTy): 
the steady state depends on the initial condition. 



FIG. 15: (a) Phase diagram at a = 10~^, for I = 0.0, ^ — 0.02, r = l/^ = 0.0, shown with the corresponding data for I = 0.008, 
^ — 0.02, r — 0.4 for comparison. Tie hnes are only shown near the critical point because for larger values of 0, there is no 
uniquely selected, smoothly banded state. This is shown in Figs, b and c. In Fig.b the steady state profiles from left to right 
at fixed ordinate correspond to tie lines left to right in the upper Fig. (a). Fig (c) shows the steady state profile in 7 (upper 
two curves) and in (f) (lower two curves) for (p = 0.16 and 7 = 4.66 with initial condition (t>{y) = + 0.4cos{ny) (solid lines) 
and with (fjiv) = 4' + 0.7cos(7rj/) (dashed lines): the "selected" state depends upon the initial condition - i.e. there is no state 
selection for I = for stresses far enough above the critical point. 



to right in Fig. |15(b) are progressively further above the 
critical point. The tie lines correponding to the smooth 
profiles near the critical point are shown in Fig. 15(a) 
alongside the correponding results at r = 0.4 for compar 
ison. Consistent with the discussion of non-universality 
in the previous section, the phase diagram for r = 0.0 is 
slightly different from that for r — 0.4 (and is different 
again from the case r = 00; not shown). 

For the spiky profiles, further from the critical point, 
the binodal o f the a ssociated tie lines is irregular (not 
shown in Fig. 15(a) I, suggesting that the steady state is 



not uniquely selected. In view of this, a natural question 
is whether selection could occur in principle (but is in- 
accessible with any realistic mesh due to the pronounced 
non-monotonicity in Wxy{y)), or whether selection can- 



not occur, even in principle. In Fig. 15(c)| we show that 
the steady state depends on the initial condition; so state 
selection appears to be lost when / = 0. This numerical 
observation is backed up by the following analytical ar- 
gument. 

In steady state, the system must obey: 
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• The force-balance equation, 

^(7, = G{(l))W^y[iTi(t>)] + 77(^)7 = S = const. (6.3) 

• The (now local) constitutive equation, equa- 
tion (EH), 

Wa(3 = Wa/iijTicf))] for aft = XX , xy ,yy . (6.4) 

• The steady-state of eauation l3.8l For the purposes 
of this analytical argument we use a simplified ver- 
sion of this equation, which we believe still captures 
the essential physics: 



= dl {A - dl^ 



(6.5) 



with 



A = m-^MMlM. (6.6) 

4> 

Integrating Ean. l6.5l twice. and using the boundary 
conditions dy4> = 0, dyCf) = for y = 0, L, we obtain 



/i = const. = A — dycf), 
where /i is an integration constant. 



(6.7) 



We now show that a solution satisfying Eons. and 16.41 
cannot in general simultaneously satisfy Eqn. 16.71 

Consider firstly Eqns. 16.31 and 16.41 Sustituting Wxy 
from Eqn. 16.41 into Eqn. 16.31 we obtain an expression for 
S'(7,(/>): this is just the family of homogeneous consti- 
tutive curves, as plotted in Fig. ^ above. Because the 
constitutive equation is local, the solution at all points 
across the rheometer cell must lie on one of these intrin- 
sic constitutive curves. Indeed, as the shear rate changes 
across the interface, the system must pass through con- 
stitutive curves of differing concentrations to maintain a 
uniform stress S. In other words, a relation <j) — 0(7, S) 
must be obeyed. The family of these curves is shown as 
dotted lines in Fig.ll6h,. For the range of stresses at which 
0(7, S) is non-monotonic, (j) must have the form shown in 
Fig. 117b in which the derivative dy(j) changes sign three 
times across the interface, as in Fig. int. [Actually, the 
forms of Fig.ll7b.c assume that the profile in 7 increases 
monotonically through the interface fFig. 117b). However 
this monotonicity will emerge self consistently from our 
argument below.] 

However we know from Eqn. 16.71 that dy4> = A — /i. A 



is plotted in Fig. I16b.c using Eqn. 16.61 together with the 
constraint 4> = 0(7, S) (imposed from Eqns. 1^31 and 
as discussed above). From this plot we see that, for any 
12, a solution that starts and ends in homogeneous phases 
(for which dycj) — A — fi — 0) 'J?] can only involve at most 
one sign change of dy(f> between the boundaries. This 
inconsistency with Fig. 117b means that a steady banded 
solution cannot exist for these stress values for which 
(^(7, E) is non-monotonic. To summarize: for stresses 



far enough above the critical point that (f>{j, S) is non- 
monotonic, a steady state solution cannot simultaneously 
satisfy Eq ns . 1^751 and 16 . 4l f which imply three sign changes 
of dyCf)) at the same time as Fan. 16.71 (which only allows 
one sign change). Therefore there a steady, smoothly 
banded profile cannot exist for such stresses. 

This argument is consistent with the sharp numerical 
which are replotted in Fig. 116b 
solution should have followed a lo- 



15(b) 



profiles of Fig. 
(solid lines): each 
cal (dotted) curve (/)(7), but instead has jumped across 
the region in which this curve is non-monotonic. (The 
stresses used to generate the homogeneous solutions 
0(7, S) in Fig. 116b were slightly different from those of 
the numerical profiles: however the trend is still clear.) 
Although Eqn. (|6.5(l is highly oversimplified, we believe 
that the failure to negotiate the interface due to the con- 
flict described above is the reason for non-selection in the 
full, numerically solved model. 

We return finally to justify our assumption that the 
shear rate must increase monotonically through the in- 
terface, and to discuss in more detail the nature of the 
banded solution when it can exist {i.e. for stresses near 
crticial point where 0(7) is monotonic). Multiplying 
Eqn. IHn across by d4>/dy, integrating on (j), and impos- 
ing = at each boundary, we find, for the simplified 
model of Eqn. 1 6'. ,51 



dcj) A]^ 0, 



(6.8) 



which is an "equal areas" construction. (<f>i and <f>r de- 
note the boundary values at ?/ = 0, L.) If, in addition, we 
were to impose that dyCf) = at each boundary, then the 
construction must automatically be as shown in Fig. 118b . 
However we did not actually impose this condition in our 
numerics, so the construction of Fig. 118b is also possible. 
This in fact corresponds to a finite system, where the 
true homogeneous state dy(j) = OVn is not quite reached 
at the boundaries. Any other equal areas construction 
(Fig. 1181 ) is not possible, for the following reason. Con- 
sider starting at point C with dy(j) — (which we do 
impose at the boundary in our numerics). Eqn. I6.7l then 
tells us that dyCj) < at this point, so the function (/)(j/) 
must curve downwards from its starting point of zero 
slope. Therefore </> locally decreases, and the system 
moves to point C . Repeating this argument, we find that 
the system can never cross to the point D. By similar rea- 
soning, the shear rate must rise monotonically through 
the interface since any initial fall (from the side of the 
low shear band) would be similarly unstable to point C 
in Fig. 118b above. 

Of course the concentration equation (|6.5() is highly 
oversimplified. For instance, a more realistic model (such 
as the one of Eqn. I3.8|l would have dependent prefac- 
tors to the dyCj) term. The equal areas result of Eqn. I\6. ^) 
is therefore specific to our oversimplified Eqn. i|6'.5|] . and 
does not hold in general. Nonetheless we believe that 
Eqn. 16.51 correctly predicts the absence of a uniquely 
banded solution for stresses far above the critical point. 
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FIG. 16; a) Dotted lines: relation between <j) and 7 for the case of a local constitutive equation, for several values of the shear 
stress, E, calculated using Eqns. 16.31 and 116.411 . Solid lines: the results of our numerics, showing that the (j>— profile cannot 
properly negotiate the interface, as described in the main text. (Note that the stresses used to generate the dotted and the 
solid lines differ slightly, but the overall trend is still clear.) b) The function A of Egn. 16.61 plotted vs. 7 using the relation of 
Fig. a), c) A replotted vs. cj>. 



via the basic conflict between the number of sign changes 
of dy<j) across the interface, described above. 

VII. CONCLUSIONS 

In this paper, we have studied the role of concentration 
coupling in the shear banding of complex fluids using the 
two-fluid, non-local Johnson-Segalman model. We have 
calculated phase diagrams for different degrees of cou- 
pling between concentration and mechanical degrees of 
freedom (molecular strain), and found a phase diagram 
qualitatively consistent with experiments on micellar so- 
lutions at dilutions well below the equilibrium isotropic- 
to-nematic transition Specific points to note are as 
follows. 

1. The coexistence plateau in the steady-state flow 
curve slopes upward with shear rate, because of 
the concentration difference between the coexisting 
bands. The overall plateau height and width de- 
crease with average concentration, terminating in 
a non-equilibrium critical point. CPCl/NaSal in 
brine [9| shows the same trend. 

2. Of the two coexisting bands, the high shear band 
has a smaller concentration due to the fact that 
concentration tends to move up gradients in the 
normal micellar strain component Wyy (where y is 
the flow-gradient direction). {W_ describes defor- 
mation relative to the unit tensor 5_, and Wyy is 
more negative in the high-shear phase than in the 
lower shear phase.) Tie lines of the phase diagram 
in the 7, 4> plane therefore have negative slope. 



3. The concentration gap is smaller for smaller values 
of concentration-coupling a oc G' {(f))/ f" {(p), and 
tends to zero in the limit a — > 0. Accordingly, the 
coexistence region of the steady-state flow curve 
becomes flat in this limit. 

4. We have described the way in which the flow phase 
diagram can be reconstructed from the family of 
flow curves E(7, (/>), measured for several average 
concentrations 4> (Fig. EU- 

5. The phase diagram and flow curves depend slightly 
on the relative size of the interfacial term in the 
viscoelastic constitutive equation to that in the 
equation that specifles the concentration dynam- 
ics. This is a concrete demonstration of how stress 
selection and the coexistence conditions of driven 
systems depend on the nature of the interface, in 
contrast to equilibrium coexistence. 

6. We find no unique state selection when there are 
no gradient terms in the viscoelastic constitutive 
equation, except for stresses that are close to the 
critical point. This implies that, for a model to re- 
produce a uniquely selected stress, it is not enough 
to simply have gradient terms only in, for example, 
the concentration dynamics. The dynamical equa- 
tions of motion for each degree of freedom must 
possess inhomogeneous terms to attain selection in 
all situations. Conversely, in situations where such 
terms are physically absent, one can expect, under 
certain conditions, no selection and hence a range 
of control parameters (shear stress or strain rate) 
for which the steady states are intrinsically history- 
dependent. 
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FIG. 17: Assuming that the shear rate varies monotonically 
across the interface (a), then for a relation 0(7, E) of Fig. ll6h , 
that is non-monotonic, the concentration must vary as in 
b), with three sign changes in 0" = d^<f) as in c). 

7. The interface width diverges at the critical point 
as a power law (S — Ec)~" with n « 0.5, although 
n differs slightly across the different order parame- 
ters. 

Although our d-JS-(/) model is highly oversimplified, we 
believe that it contains the basic ingredients required for 

I 



a first description of wormlike micellar surfactant solu- 
tions at concentrations well below the isotropic-nematic 
(I-N) transition. In particular, it incorporates the min- 
imal set of realistic degrees of freedom (tensorial order 
parameter for the micellar strain together with concen- 
tration), and unifies a non-monotonic flow curve with 
the Helfand-Fredrickson coupling between concentration 
and flow. Similar techniques could be applied to more 
involved Gates non-linear theory for wormlike micelles 

We recall a previous calculations by Olmsted et al. was 
aimed at systems of rigid rods near the I-N transition . 
In future work we hope to unify these two approaches into 
a description of wormlike micelles that is valid over the 
entire concentratio n range. This should provide a first 




FIG. 18: Of these three proposed constructions specifying 
the banded state (when it is selected, near the critical point), 
only a) and b) are consistent with the boundary conditions 
dy<l> = 0. Of these, b) is for a finite system for which dycj> 7^ 
at the boundary while a) is for the realistic physical limit 
in which the interface is narrow compared with the gap size, 
connecting two homogeneous phases in which dy(j} = OVn. 

step towards understanding the crossover regime in the 
data of Fig. |2| in which the coexistence plateau stress is 
a non monotonic function of the micellar concentration. 
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APPENDIX A: D-JS-(j!. EQUATIONS IN CARTESIAN COORDINATES 

In this appendix, we give the components of the d-JS-i/) model's equations for planar shear flow along the x direction, 
allowing gradients only in the flow-gradient direction, y, as described in Sec. IIII CI above. The x component of force- 
balance is (in the zero-Reynolds limit considered in this paper) 

= dy[G {4>) Wxy] + ■nmdy[(j) dyV^x] + Vs [ ( 1 " </>) dyVsx] ■ (Al) 

The y component of force-balance is fixed by incompressibility, V.w — 0, along with the boundary condition Vy = 0: 



= (pV,ny + (1 - (l>)Vsy 



(A2) 
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The relative velocity between the micelles and solvent (again ignoring inertial terms) is 



C 



^dy [G (0) Wyy] + 2 ^ 7? „ 9y [ dyV,ny] - 2 J^V S Oy [ (l " </>) ByV^y] - 9y | (A3) 



(1 



C 



^dy [G (0) Wxy] + ^77m 9y [ dyV„ 



1 



1 - 



:f]s dy [{I - Cj)) dyVs 



W - gd^yCb + i G' (cj)) [Wyy + W,^ ^ 111 [Wyy W + Wyy + W + 1 - W ^y")] 



(A4) 
(A5) 



The evolution of the micellar strain tensor is given by 



dtWxy + Vmy dyWxy = |(a - l)Wxx dyVmx + ^{l + a) WyydyVmx + a 



Wxy l{<Prd^Wxy 

t{4>) t{4>) 



dtWyy + 9yWyy = (a - 1) WxydyVjnX + 2 fl 

^^yydy^my ~t~ 2 3yVfi 



Wyy , l{C^fdlWyy 



t{4>) 



(0) ' 



Wxx . i{4>fdlw. 



dtWx dyWxx = (1 + a) WxydyV„ 



Finally, the concentration dynamics are 



r (0) r (0) 



= -^j^ I ^ '^'^ [G {4>) Wyy] +2rinidy[4i dyVrny] " ^ [ ( ^ " ^^^^sy] - dy^ 

I 



(A6) 
(A7) 

(A8) 



(A9) 



APPENDIX B: STATIONARY HOMOGENEOUS 
SOLUTIONS OF THE D-JS-(^ MODEL 

In planar shear, the stationary homogeneous solutions 
to Eqns. H3.7I3.9() for given 7 and are^;^.;,! = v^~v^ — 
and 



where b — l — a^. The steady state shear stress is given 



by 



Wx 



l + fe72r2(0)' 

a- 1, 



1 + a 



-Wxx = - 



1 



672 



(1 + a) 1 + 672 



W,, = Wx. = Wy, = 0, 



(Bla) 

(Bib) 
(Blc) 



T^xy = G{(t>)Wxy + (l>rim + (1 - 0)f7s7 = constant. (B2) 
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